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ABSTRACT 

The present article reviews two of the more common procedures for analyzing 
the stability and forced response of equations with periodic coefficients, 
namely, the use of Floquet methods, and the use of multiblade coordinate 
and harmonic balance methods. The analysis procedures of these periodic 
coefficient systems are compared with those of the more familiar constant 
coefficient systems. 


INTRODUCTION 

In dynamic analyses of rotating wind turbine 
systems, one frequently encounters equations of 
motion with periodic coefficients. Unlike syst^s 
with constant coefficients whose analysis techni- 
"weTT known and familiar, the analysis of 
these periodic coefficient equations are somewhat 
less familiar. The present paper reviews two of 
the more conwnon procedures for analyzing the sta- 
l)TTfty and response of these periodic coefficient 
equations, namely, the use of Floquet methods and 
the use of multiblade coordinate and harmonic 
balance methods." To put things in proper perspec- 
tive and to make comparisons, the paper will 
briefly review the constant coefficient systems 
-first”” The paper is essentially based on Appen- 
dices A, "B, C, and D of a recent report by the 
authors, (re f. 1). 

CONSTANT COEFFICIENT SYSTEMS 

Given a system of N linear differential equations 
with constant coefficients, 

M + B i + K ^ = F (t) (1) 

where M, B, and K are square matrices of order 
NxN, while £ and £ (t) are column matrices of 
order Nxl . These can be rearranged as, 
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Then, multiplying through by the inverse of the 
mass matrix gives 2N first order equations, 

y - A / = G (3) 

where A is a square matrix of order 2Nx2N, while y 
and £ are column matrices of order 2Nxl given by 
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The above rearrangement, Eq. (3), is valid provid- 
ing the mass M is not singular, which is usually 
the case with physical systems. 

(a) Stability 

To investigate stability, one sets L = 0 (which 


g^ves £ = 0) in Eq. (3), to obtain a set of homo- 
geneous equations. Then one seeks exponential 
solutions of the form y = y. Placing these 

into Eq. (3) leads to the standard eigenvalue 
probl em, 


A y = p y 


( 5 ) 


Eigenvalues pj^ of the matrix A can be obtained by 
standard numerical eigenvalue routines. If any 
eigenvalue pk is positive real or has a positive 
real part, the system represented by Eq. (3) or 
equivalently by Eq. (1) is unstable. 

(b) Forced Response 

Under steady-state conditions, the forces £ (t) on 
a rotating system tend to occur periodically in 
multiples of the rotation frequency Q. One can 
then express the force for a particular frequency 
coji, = In the form, 
io) t 

F(t) = Re(F e ) = F|^cosUj|jt - FjSlnoi^^t (6) 

The response £(t) is similarly of the form, 

1w t 

£(t) = Re(a e "' ) = q^cosw^t - ^slnco^^t (7) 

Placing Eqs. (6) and (7) into the basic Eq. (1) 
and matching sine and cosine terms gives a set of 
2Nx2N real equations. 
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where one has the matrix elements, 


G = K - H . 


H = 


(9) 


Given the amount of the m^h harmonic force present 
Fo''")and Fil"'). Eq. (8) can be solved by simple 
Tnversion* to find the response and 

for each harmonic. Then, one may sum up all 
the harmonics to give the total periodic response 
as. 


N 


(m). 


N 


a(t) = l q^-'coso) t - I qj^"’^sino)^t (10) 
m=0— m=0 — 


Finding the response £(t) this way rather than by 
direct numerical integration, allows one to assess 
the effects of a particular harmonic on the 
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resulting response of the system. 


s^tability exponent p^^ are given as 


FLOQUET METHODS 

Assume the coefficients M, K in Eq. (1) or 
equivalently the coefficients J in Eq, (3) vary 
periodically in time, rather than being constants. 
For illustrating Floquet methods, it will be con- 
venient to use the first order representation, 
namely 2N equations of the form, 

I - A(t)i^ = G(t) (11) 

where A(t) and £{t) are periodic over an 
interval t. 

(a) Stability 


The Floquet stability analysis described here 
follows that given by Peters and Hohenemser 
(ref. 2). To investigate stability, one sets 
G=^0 in Eq. (11) to obtain homogeneous equations. 
The Floquet theorem states the solution of Eq. 
(11) with £=0 is of the form 


lit) = B(t) 






where y_(t) and 
column matrices, and ^(t) 
matrix periodic over period T, that is, ^(T)=^(0). 
From the above, one can express 


{C. e } are 2N X 1 
is^a 2N X 2N square 


( 12 ) 


liO) = B(0) [C^] (13) 

Z(T) = B(T)|C|^e'’''^)= (14) 

Also, one can express y^(T) as. 


y^CO)] (15) 

y2(0)l 


where is the solution at t=T of Eq. (11) 
with £=0, for the initial conditions y] = l and all 
remaining yj(0)=0, y_(2) is the solution for y 2 ( 0 )=l 
and all remaining yi(0)=0, etc. The square 
matrix [Q] is called the "Transition Matrix." 
Equating Eq. (15) to (14) and introducing Eq. (13) 
gives. 



[Q] 


|b(0)J^c^ + |b(0)J 
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(16) 


= yJtn|X,^| = ^iln[(Ak)^ + (Ak)^] 

R I 

“k ° T ^ (20) 

I R 

The real part is a measure of the growth or 
decay of the response, as can be seen from Eq.(12). 
Values of a|^>0 (or equivalently | Xj^[>l ) indicate 
instability. The imaginary part represents the 
frequency. However, because tan“‘ is multivalued, 
one can only obtain wb to within a multiple of 27 t, 
To obtain the actual frequency and motion corres- 
ponding to the root, Pj^, one sets Ck=l and 
all other remaining Ci=0 in Eqs. (12) and (13), 
Then, using the k^" eigenvector {B(0 ))l from Eq. 
(17) as an initial condition, one would solve 
Eq. (11) with G=0 by numerical techniques for the 
resultant motion. 

Summarizing: To check for stability of a system 

of linear equations with periodic coefficients, 
obtain the eigenvalues of the "Transition 
Matrix" [Q] • If |X|^|>1, one has instability. 

The traditional stability exponent pj^ is related 
to through Eqs. (18) to (20). Two remarks on 
the above procedure should be noted. (1) The 
"Transition Matrix" [Q] can be formed by solving 
either the first order equations, Eqs. (11) with 
G=0, or the second order equations, Eqs, (1) with 
Fj=0 and periodic coefficients, whichever is more 
convenient for the integration scheme. (2) The 
above procedure will still apply even if the 
equations have constant coefficients. However, 
for such cases it is usually easier to form the 
matrix A given by Eq. (4) and obtain its eigen- 
values Pk rather than to form the "Transition 
Matrix" [Q] and obtain its eigenvalues Aj^. 

(b) Forced Response 

Solutions of Eq, (11), or equivalently Eq, (1) 
with periodic coefficients, can be obtained by 
direct numerical integration using some convenient 
integration scheme. By proper choice of the 
initial conditions, one can eliminate all tran- 
sients from the response and obtain the desired 
steady-state dynamic response by integrating 
through only one period T, instead of the very 
large number usually required to reach steady- 
state for lightly damped systems. A procedure 
for finding the proper initial conditions is 
given below. 

Solutions of Eq. (11) are of the general form, 

Z(t) = + i[p(t) (21) 


Since Cj^ are independent, one must have 

[Q] (B(0)}j^ = 

P|^t 

where A|^=e are the eigenvalues of the [Q] 
matrix. One then has the relation 

P[^ ~ j" ^ ^ ^ ^k 
from which the real and Imaginary parts of the 


where homogeneous solution and yp(t) 

is the particular solution. One can obtain a 
complete solution of Eq. (11) numerically for any 
given set of initial conditions. Call this solution 
y^(t). One can add any number of additional homo- 
geneous solutions Ay^(t) having different Initial 
conditions, to this solution. This would give a 
new solution to Eq, (11), 

il(t) = (22) 
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which would have different Initial conditions 
than those for ^^(t). 


One can obtain all the homogeneous solutions of 
Eq* (11) by solving Eq, (11) with G=0 a total of 
2N times, subject to the initial conditions y]=l 
and all remaining y-j=0, then y 2 =l and all re- 
maining y-}=0, etc. In fact, this was done 
earlier to Investigate stability and resulted in 
the 2N homogeneous solutions 
etc., respectively. ^ \ L ^ 

Thus one may write 






(23) 


[Q(t)] 


Where [Q(t)j is the transition matrix at any in- 
stant of time, and are 2N arbitrary 
constants.. The new solution Eq, (22) can be re- 
written as 


^ W = ^£(t) + [Q(t)]C (24) 

For a periodic solution over period J=Z-n/U, one 
miist have ^(T)=^(0). Placing Eq. (24) into this 
condition and solving for the arbitrary constants 
C gives, 

Z£(T) + [Q(T)]C = i^{0) + [Q(0)]C 

C = [l - [Q]]-’ {ze(T) - yc(0)} (25) 

where it was noted that [Q(0)]=1_, and [Q(T)]=[Q] 
is the "Transition Ffetrix” found earlier for the 
stability Investigation. Placing these values of 
C back into Eq. (24), the initial conditions for 
insuring a periodic solution become 

LW = yp(0) + (i^£(T) - ^£(0)} (26) 

One can then solve the basic Eq. (11) numerically 
with these initial conditions to obtain a 
periodic solution over one period. It should be 
noted that If one had chosen the Initial con- 
ditions for yr(t) as yr(0)=0, one would obtain 
simply 


x(0) = ygd) (27) 

This is a particularly convenient form for finding 
the initial conditions for periodic solutions. 

An alternative form for determining the proper 
initial conditions for periodic solutions has 
been proposed by Friedmann and his coworkers 
(refs. 3 and 4) In their work on wind turbines, 
namely. 


^(0) = [i-Q]”' a / [Q(t)]'^F(t)dt (28) 

0 

This Is similar to Eq, (27), but does not use 
It seems easier to obtain yr(T) with initial 
conditions y^(0)=0 and use ^ Eq. (27), rather 


than obtaining [Q(t)] at every point and perform- 
ing the indicated operations required by Eq,(28), 

The general procedure described by Eqs, (21) to 
(27) may be extended to deal also with nonlinear 
equations. 


i - A(t)y^ = F_(t, y_) (29) 

where the right hand side now contains nonlinear 
functions of the coordinates. An iterative 
variation of the previous linear procedure to ob- 
tain the initial conditions for periodic solutions 
of nonlinear equations is used by the MOSTAS Code 
(ref, 5). The procedure is as follows. First, a 
numerical solution yr(t) is obtained to the non- 
linear Eq. (29) for ^some estimate of the initial 
conditions yr(0)* Then each of the 2N elements of 
yp(0) is ^ perturbed a small amount and the 
resulting 2N solutions are obtained. This in- 
volves solving the nonlinear Eq. (29) subject to 
the Initial conditions, 
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(30) 


and will result in 2N responses of the form 


z^’^(t) = )^£(t) + A^y^(t) (31) 


where A ^(t) represents the effect of each per- 
turbation c., and is found by subtracting yp(t) 
from each ^ of the 2N resulting responses 
y^(^)(t). One can then express the total solution 
approximately as. 


lit ) . i^it) * 
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(32) 

which is in the same form as Eq. (23). Then, 
again requiring the periodicity condition y(T)= 
y(0) and following through as before, will result 
in the same relation Eq, (26) found previously^ 
Because of the nonlinearities now present, the 
elements of [Q] as found from Eqs. (32), (31), (30) 
may vary with the amplitude of the initial con- 
dition used, yfr(0)+€. . This is in contrast to 
the linear case where [Q] remains always constant. 
Hence, an iterative application of Eq. (26) with a 
new corrected £r(0) should be done* If the non- 
linearities arF^not too great, convergence to the 
required ^^(0) should be rapid. 


It should be remarked that the numerical procedure 
for forced response described in this section, can 
also be used for the constant coefficient linear 
case, although it is probably easier there to ob- 
tain the solution by using Harmonic response 
methods given by Eqs. (6) to (10). However, for 
cases where there is some nonlinearity, the 
present iterative approach becomes attractive. 
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MULTIBLADE COORDINATES AND HARMONIC BALANCE 

Given a rotor with N blades rotating with rotation 
speed n, attached to a flexible tower. Because 
the tower motions x. are described in a fixed 
reference frame while the blade motions 3- are 
described relative to a rotating frame, ^the 
resulting equations may have mass, damping, or 
stiffness coefficients which are functions of the 
azimuthal position of the blade ip. . A typical 
such set of equations is given, for example, in 
refs. 6 and 7 as, 

2 k 

Mx + C X + k X + S 4). = F (t) 

X X dt"^ K X 

Sx cos 4;,^ + 

(k = 1, 2, ... N) 

where the azimuthal position 'i^|^ is, 

4 ),^ = nt + (k - 1) 2 tt/N (34) 


Mx + *^x^ f ^ '’c " ° 

I [l bs + Cgb^ (kg - I«2)b3 

-2iH6^-£lCg bj = 

j jsx + 2 m 6^ + ncgb^ + li)’^ 

* ‘A * <*^6 - 

" ['^o * S'o * Vo] ■ " 

These equations may then be investigated for 
stability using the standard constant coefficient 
techniques described earlier. For additional de- 
tails and applications of multiblade coordinates, 
see Hohenemser and Yin (ref. 8). Multiblade co- 
ordinates were originally Introduced by Coleman 
and Feingold (ref. 9) in their studies of heli- 
copter ground resonance. 


0 

(37) 


= 0 


The first equation above represents force equi- 
librium for the tower motion x, while the remain- 
ing N equations represent force equilibrium for 
the motion of each of the N blades , The 
above equations are readily generalized to more 
tower motions x., and more blade coordinates for 
each blade 

(a) Stability 

To examine Eqs. (33) for stability, one sets 
F^=0 and F^(^)=0 to obtain homogeneous equations. 

For rotors with 3 or more blades N>3, one may 
eliminate the periodic coefficients in these 
equations by introducing new multi blade coordinates 
bo(t), bg(t), b^(t) such that 

= bj^(t) + b^(t) sin 4 ^ 1 ^ + b^(t) cos i|;|^ (35) 

Substituting these into Eqs. p3), then multiply- 
ing the last N equations by sin , cos , and 1 
respectively, then summing these ^last N ^equations 
and noting that 

N N 

I sin = I cos ijj. = 0 
k=l k=l ^ 

N 2 2 

I sin ipK - I cos"^ = N/2 
k-1 k=l ^ 

N 

I sin i|). cos = 0 

k=l ^ 

results in a new set of differential equations in 
the variables x, b , b , b which now all have 
constant coefficients , ^namely. 


For N 3 

(36) 


For rotors with 2 blades, N=2, the analysis is 
more difficult because the rotor disk no longer 
has polar symmetry. If the same multi blade co- 
ordinates given by Eq.(35) are used in the basic 
Eqs. (33), the periodic coefficients would not be 
entirely eliminated since now, 


^ 2 

I ^ 1 , = 1 - cos 2\p. 
k=l ' 


c 2 

I cos = 1 + cos 2iJj. (38) 

k=1 ^ ‘ 


2 

I sin % cos % = sin 2^;, 
k=^ K K I 

instead of the convenient constant terms given by 
Eqs. (36). A rough estimate of the stability and 
response can be obtained by simply time-averaging 
the resulting cos 2\p, and sin 2ij^. variations to 
zero and using only the constant coefficient 
terms. This is equivalent to setting N=2 in the 
multiblade transformed Eqs. (37). 

For more accurate estimates for these 2-bladed 
rotors, one may use harmonic balance methods. 

This consists of first introducing new coordinates 
b-j-(t) and b^(t) for these two blades such that, 

B^’^ = Bj + 6^ , B^^^ = 

then summing and subtracting the last two equations 
of Eqs. (33) while noting that sin 4 ; 2 "“Sln . and 
cos i|/ 2 =-cosi/).| , then expanding each '^of the coor- 
dinates in a harmonic series, 

x = Xo + x-|^s1ns^t+ x-j^cosJlt+ X25Sin2S7t+ ... 

bj = bjQ + b-p^^sin ?2t + b^^^cosl^t + .** 

‘’a ^ '^AO b^^ ^cos^lt+ ... (40) 
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wh 6 rfi Xp ) c * • Jt 9 ST'S si 1 "functions of 

time* ^ ^ ^ ^ Plscing these into 

Eqs.(33),and bslsncing out each hannonic term In 
each equation will yield a truncated series of 
constant, coefficient differential equations. 

These aquations may again be examined for stability 
using the standard constant coefficient tech- 
niques, described earlier. 


R.H.S. = 


f F^(t) 


I Sin 


k=l ^ 

N fkl 

I Fg' (t) cosi];. 
k=l 


(42) 


Often, depending on the form of Eqs.(33), the re- 
suiting constant coefficient differential equations 
will uncouple into several smaller coupled systems 
of equations which may be examined independently 
of one another. For example, for the case of 
Eqs.( 33 ), one smaller coupled system would involve 
the variables Xq, X 255 ^AlC* ^AIS’ while 

another would Involve x^^, x^^, b^Q, by^ 2 Qf ^a2S ’ 

,,, For such systems, one could use an alternate 
extended form of the multi blade coordinate trans- 
formation Eq.( 35) namely, 

X = Xq + y^2S 2^1t+ x^^cos 2Qt+ ... 

= b-j^ simple + b.j^ cos^^ t ... (41) 

together with the harmonic balance method to solve 
the problem. This works here, since the form 
given by Eq.(41) exactly duplicates the motion of 
the two blades given by the general case Eqs,(39) 
and (40), since sin i|; 2 =-sin , cos i /^2 » 

and only Xq, X 25 , X 2 Q, ^AIS’ "* would be 

present; However, in more general cases (for 
example. If the first equation of Eqs.(33) had an 
additional term M-jX cos i|^.| or k.|X cos present), 

the resulting equations would not split into two 
smaller groups, and the general harmonic balance 
method Eqs.(39) and (40) would have to be used. 

Indeed, for the more general case mentioned above, 
one would also investigate the system for direct 
Mathieu equation type instabilities of half 
Integer order Q/2, 3fi/2, ... by introducing 
additional harmonic terms sinmQt and cos mQt 
where m=l/2, 3/2, 5/2, ... into Eqs.(40), and 
harmonically balancing as before. These terms 
would not couple in with the previous equations 
and can be solved independently of them. The 
primary instability region would result from the 
Q/2 terms. See Bolotin (ref. 1C? for further de- 
tails of the general harmonic balance method. 

Also see Sheu (ref. 7) for an application of the 
alternate extended form of the multi blade trans- 
formation Eq,(41), to a simple two bladed rotor 
in ground resonance, 

(b) Forced Response 

For rotors with 3 or more blades, N^3, one uses 
the multiblade coordinate transformation Eq.(35) 
to eliminate the periodic coefficients in the 
basic equations of motion Eqs.(33), as described 
in the preceding section. The equations then re- 
duce to the constant coefficient equations given 
by Eqs.(37), only now the right-hand-sides are 


l F^^\t) 

lk=l ^ J 

instead of the previous value of zero. Under 
steady-state conditions, the tower and blade forces 
generally occur periodically in multiples of the 
rotation frequency and can generally be 
expressed as. 


= F^o + F^-i^sin + F^iccosg.^ 
+ 2 .1,, + ... 


(43) 


where i|;. = Qt + (k-l)2iT/N. Placing the above 
forces ^ Into Eqs.(42) and using the trigonometric 
identities and summations, 

sin sin 4 ^l^ = ^ cos(m-l cos(m+l 


COT m\pj^ sin^ = etc. 

N fN sin mi|j, 

I sin m^. = 

k=l ( 0 

N fN cos mp. 

I cos 1 #. = \ 
k=1 1 0 


m = N, 
■t- m N, 
in = N, 
in N. 


2N, ... 
2N, ... 
2N, ... 
2N, ... 


(44) 


one can obtain the right-hand-sides of Eqs,(37) 
in terms of either constants or harmonic functions 
of ntlt. The forced responses x(t), b (t), b (t), 
b (t) can then be found using the standard 
techniques for constant coefficient systems dis- 
cussed previously. It should be noted that be- 
cause of the multi blade transformation Eq. (35), 
the resulting responses for the tower motion and 
blade motions corresponding to the harmonic 
oijyj = would be of the form, 

X = Xn coso)^t - X, sinca^t 
R m I m 

cosco^t - bgj Sln6i^t (45, 

+ ^‘’SR 

+ (b(.R cosio^t - bj-j sinoi^t) cos.),,^ 

The tower thus oscillates at frequency (orn In the 
fixed frame whereas the blades may oscillate at 
frequencies w , - Q relative to_ the 

rotating frame. 

For rotors with 2 blades, N = 2, the multiblade 
coordinate transformation Eq. (35) does not eli- 
minate the periodic coefficients, but rather 
changes the cos variations to cos 2 ipj^ varia- 
tions. A rough estimate of the response can be 


169 



obtained by simply time-averaging the resulting 
sin and cos variations to zero, and then 
proceeding with the remaining constant coefficient 
terms, as was done for the N>3 case. The results 
are likely to be somewhat off for the second har- 
monic, sin 2\p^ and cos responses. 

For more accurate estimates for these 2-bladed 
rotors, one can use the harmonic balance methods 
of the previous section. The steady-state peri- 
odic tower and blade forces given by Eqs, (43) 
are substituted into the basic equations of 
motions Eqs, (33). One then introduces the new 
coordinates given by Eqs. (39), then sums and 
subtracts the last two blade equations, then 
expands the tower and blade motions as given by 
Eq, (40), only now the coordinates xg* xjc, bjo> 
^TIS* by\Q,.,,etc, are taken to be constants rather 
than functions of time. Harmonically balancing 
the various terms in each equation results in a 
truncated set of algebraic equations which can be 
solved to obtain the coordinates Xq, xig, 
etc., corresponding to the given forcing excita- 
tions ^ 30 * F 3 is,-..etc. The resulting 

tower and blade motions are then given directly 
by Eqs. (40) and (39). The resulting set of al- 
gebraic equations will often uncouple into smaller 
coupled sets of equations which can be examined 
independently of one another. This procedure is 
similar to that for the constant coefficient 
forced response case Eq. (8), except now, the 
periodic coefficients couple the different har- 
monics together. Thus, the solution will consist 
of many harmonics wQ, even if only one forcing 
harmonic F^^^ were present alone, 

ROTATING COORDINATES 

As an addendum to the previous multi blade coordi- 
nates and harmonic balance methods, it should be 
mentioned that for some problems, the use of ro- 
tating coordinates is also convenient. For ex- 
ample, in the case of a 2-bladed rotor on isotropic 
tower supports (same tower mass, damping, and 
stiffness In two directions, X] and X 2 ), Eqs. (33) 
would read, 

2 k 

(46) 

Me Me ■ ^ ^ ^ 

Sx^ cosif'i^ - Sx’2 

= (k=l,E) 

One can then express the tower motions x] and X2 
in terms of rotating coordinates and C2 which 
rotate with the rotor, as 

x^ = ^ cosQt + ^2 sinJ2t (47) 

X 2 = -C-j sinQt + ^2 cosCIt 

where the rotation i|;-|=f2t is taken from the ^2 


axis towards the X] axis. Placing these equations 
into Eqs. (46), then multiplying the first two 
equations by cosip-j and sinip] respectively and sub- 
tracting, then multiplying the first two equations 
by simjji and cosij;i and adding, then subtracting 
the third and fourth equations, then adding the 
third and fourth equations will result in a new 
set of differential equations in the variables , 
b^, bj which now all have constant coefficients, 
namely, 

+ En?2 - + RC 2 ) + 

+ ES(B^ - = F^^cosnt - 

M(C2 - + C^(^2 - + ^^2 

- 4SQB^ = F^^sinJ^t + F^2'^osQt 

2S{'t^ + EJ2C2 ■ ^ ^'^e^A 

_ F 0) F 

B 6 

Eli^ + ECgBy + EkgBj = FgO) + Fg(2) (48) 

In the above, Bj=(B^^ )/2 and 6«=(B^^ ^-6^^b/2 

are the same coordinates introduced earlier in 
Eqs. (33). These differential equations may then 
be investigated for stability and forced response 
using the standard constant coefficient techniques 
described earlier. Such analyses of a 2-bladed 
rotor on isotropic tower supports were also per- 
formed by Coleman and Feingold (Ref. 9) in their 
studies of helicopter ground resonance. 

Rotating coordinates are often used in rotating 
machinery shaft critical speed problems, and are 
useful for dealing with problems of rotors with 
unsymmetrica! mass, unsymmetrical damping, or 
unsymmetrical shaft stiffness supported on iso- 
tropic bearings. See for example, Bolotin (Ref, 11). 
For such problems, one can readily set up the 
equations of motion in the rotating frame direc- 
tions, and the fixed supports will introduce no 
periodic terms because of their isotropic nature. 

For vertical axis wind turbines, such rotating 
coordinates for the blades are useful since the 
tower supports are generally isotropic due to the 
symmetrically arranged guy wires. For horizontal 
axis wind turbines, the tower supports are generally 
not isotropic, hence periodic coefficients will 
remain in the equations when using rotating coor- 
dinates. If the support anisotropy is not too 
large, one can again additionally introduce har- 
monic balance methods to eliminate the periodic 
coefficients, as was done in the previous section, 

CONCLUDING REMARKS 

The present article has reviewed two of the more 
common procedures for analyzing the stability and 
forced response of rotating systems with periodic 
coefficients, namely, Floquet methods and multi - 
blade coordinate, harmonic balance methods. Also, 
the use of rotating coordinates \as discussed. 
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The Floquet methods are based on a convenient 
numerical integration scheme and involves the com- 
putation of the "Transition Matrix," [Q], from 
which stability and the initial conditions for 
steady-state response solutions can be obtained. 
These methods seem attractive for large systems 
and can be modified to include non! ineari ties in 
the equations. 

The multi blade and harmonic balance methods in- 
volve first the Introduction of multi blade coor- 
dinates in order to take out the periodic coeffi- 
cients from the blades [Eqs.{35) for N>3], or to 
obtain a better ordered system of equations 
[Eqs,(39) for N=2]. Then, harmonic balance 
methods Eqs.(40) are used to deal with any remain- 
ing periodic coefficients. These methods seem 
attractive for smaller systems and can give con- 
siderable insight into the origin and nature of 
instabilities and the various harmonics present 
in the forced response. 

Rotating coordinates can also be used to effectively 
eliminate the periodic coefficients in problems 
involving unsymmetrical rotors on isotropic tower 
supports. These can often be used in rotating 
shaft critical speed problems and for vertical 
axis wind turbines. If the support anisotropy is 
not too large, harmonic balance methods may 
additionally be used to deal with any remaining 
periodic coefficients. 
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QUESTIONS AND ANSWERS 
J. Dugundji 


From: P,R. Barnes 

Q: The Singularity Expansion Method (SEM) introduced by Dr. Carl Baum of Kirkland AFB, 

Albuquerque, NM is another, perhaps better, approach to solving these problems. 

Do you know about SEM? 

A: I do not. I have just dealt here with two of the more common methods for deal- 

ing with these problems . 

From: W.E. Holley 

Q: Are you aware of any treatments of stochastic problems with periodic coefficients? 

A; I have not dealt with that aspect of the problem ^ so I am not aware of them. I 

believe though that there is considerable literature on that subject. 
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